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Abstract 

We review the effect of uniaxial strain on the low-energy electronic dispersion and Landau 
level structure of bilayer graphene. Based on the tight-binding approach, we derive a strain- 
induced term in the low-energy Hamiltonian and show how strain affects the low-energy electronic 
band structure. Depending on the magnitude and direction of applied strain, we identify three 
regimes of qualitatively different electronic dispersions. We also show that in a weak magnetic 
field, sufficient strain results in the filling factor v — ±4 being the most stable in the quantum 
Hall effect measurement, instead oi v — ±8 in unperturbed bilayer at a weak magnetic field. To 
mention, in one of the strain regimes, the activation gap at = ±4 is, down to very low fields, 
weakly dependent on the strength of the magnetic field. 
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I. INTRODUCTION 



Electrons in monolayer graphene are chiral quasiparticles with Berry phase tt, Dirac-like 
linear dispersion and Landau quantization in magnetic field leading to a peculiar spectrum 
with a Landau level at "zero energy" (the Fermi point separating conduction and valence 
bands in a neutral structure) [1-4J. Electrons in bilayer graphene also exhibit an exotic 
behaviour. They can be viewed as 'massive chiral fermions', that is, quasiparticles with 
a parabolic spectrum but also Berry phase of 27r [1] HHB], and their high-magnetic-field 
spectrum features an eight-fold degenerate Landau level at zero energy. It has been noticed 
that the topology of the constant-energy lines in the dispersion bands in bilayer graphene 
may change, from single-connected, almost circular at higher energies, e ^ lOOmeV, into 
several disconnected at lower energies e ^ ImeV: a transformation known in the physics of 
metals as the Lifshitz transition [7J , which occurs when the Fermi level passes the saddle point 
in the spectrum. Theoretically, details of this transformation in the spectrum and the energy 
e* at which topology of the energy band changes are determined by the interplay between 
the nearest-neighbour and skew interlayer hoppings of electrons. For energies \e\ <C the 
single-particle electron spectrum consists of four Dirac cones with a linear dispersion, three 
with a Berry phase tt and one with — tt. In this article, we show that the topology of the 
low-energy single-particle dispersion of bilayer graphene critically depends on mechanical 
deformations of the crystal. Strain not only infiuences the energy e* at which the topology 
of the dispersion changes but more critically, determines the number of Dirac cones in the 
low-energy part of BLG spectrum, from four in the unperturbed bilayer, down to two (both 
with the Berry phase equal to tt) in a strongly strained crystal. To mention, no such dramatic 
changes occur in strained monolayer graphene where deformations of the crystalline lattice 
induce merely a shift in the position of the Dirac point in the momentum plane [8^ 9j. We 
also track the strain-induced spectral changes in the bilayer down to the evolution of the 
Landau levels for electrons, which enables us to predict what features in the quantum Hall 
eflFect in bilayer graphene crystals would persist down to the lowest magnetic fields and 
lowest carrier densities. 
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FIG. 1: (a) Side view of the unperturbed bilayer graphene lattice. The top (bottom) layer is shown 
in grey (black), (b) Electronic band structure of unperturbed bilayer graphene in the vicinity of 
the Brillouin zone corners K and K\ In the top left is shown a detailed view of the low-energy 
band structure around the Lifshitz transition at energy e*. (c) Side view of the deformed bilayer 
graphene lattice. The lattice was deformed using parameters S = 0.20, 6' = —0.04, 9 = 25° and 
5r = ^(l,l). 

II. TIGHT-BINDING MODEL FOR STRAINED BILAYER GRAPHENE 



Bilayer graphene consists of two honeycomb layers of carbon atoms arranged according to 
Bernal stacking [lOj (see Fig. [l](a)) with four inequivalent carbon sites in the unit cell, Ai^ Bi 
in the bottom layer and ^42, B2 in the top layer, and with A2 positioned over Bi. Hexagonal 
symmetry in the plane leads to three nearest neighbours at the same distance surrounding 
each site, their position described with the help of vectors ei = r^B(0, 1), 62 = — |) 

and 63 = VAsi^^ — |), where Tab is the carbon-carbon distance. 

For the unperturbed bilayer, the electronic spectrum is gapless, with the conduction and 
valence bands touching in the corners of the hexagonal Brillouin zone, denoted as the K 
points at the energy conventionally used as the e = reference level, which also coincides 
with the Fermi level in undoped bilayer graphene. Two additional bands are split from the 
neutrality point by energy e ^ ±71 (see Fig. [l](b)). At very low energies, e ^ ImeV, the 
isoenergetic line undergoes a topological transition as it splits from a single-connected line 
into four disconnected parts, each resulting from a separate Dirac cone, as shown in the top 
of Fig. |3b) . 

In the suspended graphene devices pTHTS] deformations can develop in the atomically 
thin carbon membranes due to the motion of contacts as the sample is cooled down. Un- 
derstanding the influence of deformations of bilayer graphene in particular on its electronic 
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properties is, now, becoming a pressing matter, in order to distinguish strain effects in sus- 
pended devices from the spectral changes expected in the recently predicted low-temperature 
phase transitions induced by the electron-electron interactions in undoped bilayer p!6l - [2T] 
which are being searched for in suspended bilayer graphene devices [ISl-fTS]. Below, we char- 
acterise strain using angle 6 between the principal axis of strain tensor and coordinates in 
Fig. [H and its eigenvalues, 5 and 5' . Recently investigated interlayer shear deformations 
[22j can be included by allowing a shift of the top layer with respect to the bottom one, by 
Sr = {Sx^Sy). As in monolayers, stretching of the lattice changes its symmetry and makes 
the intralayer hops Al-Bl {A2-B2) dependent on the direction of the hop, from a single 
value, 7o in the Slonczewski-Weiss-McClure parametrisation f23l [24]. to three different val- 
ues for the hops along bonds ei, 62 and 63, 70^^ = 7o 1 + ^0 (^T^^^'^ + ' where 
n = 1, 2 or 3, and I = (sin 20, cos 20) takes into account unilateral deformation while the 
second term, ^i^, accounts for the 'hydrostatic' rescaling of the lattice period and leads 
only to a change in the value of Dirac velocity, v ^ v [l + r]Q^^]. Here, tjq = ^^^^ 
quantifies the change of the intralayer A-B hopping upon the change of the distance of car- 
bon atoms on the plane. To describe the effect of strain in bilayers, we need to account for 
additional changes in the interlayer hoppings. The direct A2-B1 interlayer coupling, 71, may 
be changed by shear, 71 ^ 71 + 0((^^), which has no bearing on the topology of electron 
bands at low energies. More importantly, strain and shear make the value of the skew in- 
terlayer coupling 73 dependent on the direction of the hop, leading, again, to three distinct 
couplings 7j") = 73 {1 + [f^'{^l - ^) + ^] ^3}, where 773 = for the three 
hops indicated in Fig. [TJwith black dashed lines. 

The variation of the hopping elements in the tight-binding model by strain and the 
lower symmetry of the bilayer lattice result in modification of the dispersion relation for 
electrons at low energies. The latter is most convenient to capture using an effective 2x2 
Hamiltonian, which describes in details the pair of low-energy bands near the Dirac point in 
bilayer graphene corresponding to electron states located predominantly on the sublattices 
Ai and B2 (those coupled by the skew hop 73, but not by the direct hop 71). The calculation 
leading to such effective Hamiltonian includes two steps. One consists in the Schrieffer- Wolff 
transformation [25j projecting the four-band model onto the effective two-band model [5j. 
The other takes into account that cancellation between the contributions of the 70 and 73 
hops in directions ei, 62 and 63 for the coupling of plane wave Bloch states on ^41/2 and B1/2 
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sublattices in the corners K and of the BriUouin zone is no more exact. Note that lattice 
deformations also distort the hexagonal shape of the BriUouin zone, what is automatically 
taken into account in our analysis. As a result, the hopping integrals between sublattices 
are modified from ^v{px — i^Py) to ^v{px — i^Py) + ^he intralayer A-B hops (where 

V = and ^0 = f (^^ - ^^Oe"^'%o^o), and from ^^3^^ + ^^Py) to ^^3^^ + i^Py) + A3 for 

the interlayer skew hops (where = ^^^^ and A3 = \{5 — (^')e~^^^737y3 — f^i^e^'^Ta^a with 
Lp denoting the anticlockwise angle between negative direction of the y axis and vector 5r). 
Then, the resulting two-band model can be described using the Hamiltonian 

71 y (e^TT + ^0)' I \ e^37rt + ^* 

where = ^ and tt = Px + Wy Following the suggestion [8j that in monolayers the 
effect of homogeneous strain is equivalent to a constant vector potential, which can be 
eliminated by a gauge transformation equivalent to a shift of the valley centre from the 
BriUouin zone corners [K (^ = +1) or (^ = —1)], we introduce here valley momentum 
p = p+(^^^^^(cos 20, — sin 20), and arrive at the effective Hamiltonian for strained bilayer 
graphene: 




where 



3 , X _,-9/Q / ^ 3 • 5r 

2 Tab 



^ = 7(% - ^0)736 ''^^S - S') - -73%e'^ — , n =px + ipy 



III. THE INFLUENCE OF THE ELECTRON-ELECTRON INTERACTION ON 
THE STRAIN-INDUCED TERM IN THE LOW-ENERGY HAMILTONIAN 

In this section, we investigate the influence of electron-electron interactions on the low- 
energy Hamiltonian of (strained) bilayer graphene using the renormalisation group (RG) 
approach. The derivation of the RG equations is based on the analysis of the leading 
logarithmic divergences in the perturbation theory and their subsequent resummation. It 
follows the line of Ref . [19] , though in this paper we will not consider the renormalizations 
of the marginal short range interaction giving rise to the instability towards nematic state 



[H]. As usual, the bookkeeping of the perturbation theory is greatly simplified by using the 
diagrammatic technique. 

We pictorially represent the Hamiltonian ([T]) as —H = Q + A + \I\; 
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= i ° '^"^ I; A = -^.3| Tt |;n = -| " " I- (2) 



In the energy region where the RG fiow occurs, two last terms in Eq. ^ can be considered 
as the perturbation so the Green function has the form 

- =G{e,p) = -^. (3) 

The Coulomb interaction denoted by 

'v^N^N^N^ _ 27re\ 

^,9 ~ g ' 



(4) 



is divergent at g ^ so all the terms in 1 /N approximation should be summed which gives 

— ■ — o- • " 

The explicit calculation of the polarization operator yields 



, , 4a;^ + 4\ 2arctana; — arctan(2a;) 
In . o . - + 



4a;2 + 1 y X 



(6) 

where = 4 is the total degeneracy of the electronic spectrum. All the calculation is 
controlled in 1/A approximation. Using Eqs. ^ and ([g]), we obtain 



^^^^^ = ^ ^ { ) 



I 2m(jO \ 



+ n Nm 



(7) 



The last equation is written under assumption that ga^ <C 1 ( = 1/Nme^ being the 
screening radius), which is always true in the region of applicability of the Hamiltonian (1) 
of the main text. 

To derive the RG equations we, first, calculate the correction to the self-energy. On each 
step of the RG the integration over uj along the interaction line is restricted by — ^E" < 
|cc;| < E. Its energy derivative gives the correction to the quasiparticle weight 

de (8) 



which can be always ehminated by the rescahng of the Fermionic operators and thus does 
not affect by itself any gauge invariant quantity. We will perform the rescaling in such a way 
that the quasiparticle weight is set to unity, so that all the corrections have the meaning 
of the observable spectrum (note in passing that all the corrections to the observables are 
logarithmic whereas SZ is divergent more severely). As a bonus, the scalar vertex is also 
not renormalized due to the gauge invariance 






The correction to the effective mass is found from 
The other corrections are 

A fdE\ dx D{x){l - 3x^) 



(9) 





(10) 



Calculating the integrals in Eqs. ([9])- (10), substituting = 4, and dX = dE/E^ we 
obtain: 

= 0.11^, dxm-^ = -0.02m-\ dxVs = -0.02^3, (H) 

where A = In ^ and e is the running energy scale. The electron-electron interaction enhances 
the strain-induced term stronger than other parameters, and at energies e ^ where the 
influence of strain plays a dominant role in determining the electron spectrum, we substitute 
\w\ exp(0.11 In ^) ^ \w\^-^^^^-^^ ^ in the Hamiltonian ([T]). 



IV. LOW-ENERGY ELECTRONIC BAND STRUCTURE OF STRAINED BI- 
LAYER GRAPHENE 

The resulting low-energy electronic dispersions of the conduction band [the band structure 
described by Hamiltonian in equation ([T]) is always electron-hole symmetric] of strained 
bilayer graphene in the vicinity of the valley are displayed for representative values of 
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FIG. 2: Three-dimensional plots of the low-energy electronic dispersions in the conduction band 
of (strained) BLG, for representative points in the ^w) space. The energy and momentum 
are measured in the units of e* = mv\l2 and mvo^^ respectively. The following values were used to 
obtain the graphs: the effective mass m — O.OSSmg {m^ is the electron mass), vo^ — lO^m/s. The 
resulting energy e* ^ ImeV. 



^w) in Fig. 2, The dispersion for the case of the unperturbed bilayer, that is for it; = 0, 
is shown in Fig. 2(a). It is quasi-parabolic at energies e e* = with '^3 responsible 
for the trigonal warping of isoenergetic lines. At e = e*, the isoenergetic line undergoes a 
Lifshitz transition and splits into four Dirac cones [5j. Remaining dispersions in the top 
row depict the bottom of the conduction band for '^w ^ 0, = 0. As the value of "Rw is 
increased (Fig. [2](b)-(c)), the two side cones positioned off the axis and the central cone 
move closer, as shown in the graph for w = . Those three cones collide for w = Se*, and 
for w > 3e* only two Dirac cones remain, as shown in Fig. [2](c) for w = 5e^. If < 
is negative (graphs (d) and (e) in Fig. [2]), the central cone and the side cone positioned on 
the X axis approach each other and collide for w = — e*, creating a local, quasi-parabolic 
minimum. This minimum persists for some range of the strain, although it lifts off the 
e = plane, as shown in the graph [2](d) for w = —56*. Eventually, for ?R:w = —96*, the 
minimum merges with a saddle point and with further decrease of again only two Dirac 
cones remain in the spectrum, as in Fig. [2](e) for w = —9.56*. A contrasting situation of 
= 0, ^ 0, is presented in the graphs in Fig. [2](f)-(h) (due to symmetry, dispersions 
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FIG. 3: Diagram showing what distinctive features are present in the low-energy electronic disper- 
sion of BLG as a function of the deformation. 

for opposite values of are mirror reflections of each other with respect to the Py = Q 
plane and hence only situation of < is described here). Again, starting from the 
unperturbed system and increasing the magnitude of leads to two of the cones moving 
closer to each other (Fig. [2](f) for w = — ze*) and colliding with a creation of a local, quasi- 
parabolic minimum which lifts ofl[ the e = plane (graph (g) in Fig. [2]for t(; = — i3e*) and 
eventually disappears, leaving only two Dirac cones (graph (h) in Fig. [2] for it; = —ibe"). We 
identify three qualitatively different regimes in the {^w^^w) space: (i) spectrum contains 
four Dirac cones, (ii) spectrum contains two Dirac cones and a local minimum, created by 
a collision of two of the cones and which in general does not touch the e = plane, and (iii) 
spectrum contains only two Dirac cones. The extent of those regimes is shown in Fig. [3| 

Strain-induced deformation of the low-energy electronic dispersion gives rise to distinctive 
features in the low-energy density of states in each of the strain regimes. In particular, 
existence of a local parabolic minimum results in a sharp step in the density of states, while 
saddle points lead to van Hove singularities. Note that for t(; 7^ 0, the symmetry between 
the side cones is broken and the saddle points are no longer all found at the same energy 
e*, leading to more than one van Hove singularity and more than one Lifshitz transition. 
Densities of states representative of all the regimes are shown in Fig. [4] for w = w = — ie*. 
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FIG. 4: Low-energy density of states in the conduction band of (strained) bilayer graphene for 
representative values of w. 

w = — i3e* and w = — i56*, that is points that he on the = hne. Indeed, for it; = 
only one van Hove singularity is present, whereas three can be traced for w = — ie*. A sharp 
step, corresponding to the contribution of a quasi-parabolic part of the band structure to 
the density of states, is visible for w = — zSe*. It is followed by a van Hove singularity, as 
for this dispersion the saddle point lies not far above the position of the local minimum. 
Finally, for w = — i5e*, the sharp step is washed out (as the local minimum merged with the 
saddle point) and only one van Hove singularity remains. 



V. LANDAU LEVELS OF STRAINED BILAYER GRAPHENE 

Presence of strain in bilayer graphene may also lead to new features in the low-energy 
Landau level (LL) spectrum. To obtain Landau level spectra for bilayer graphene in a 
perpendicular external magnetic field, we use the two-band Hamiltonian ([T]), where the 
electron momentum, p = —ihV — eA, incorporates the electromagnetic vector potential 
A. For simplicity, we use the Landau gauge, A = (0,5x). The operators tt and iv fulfil 
the role of raising and lowering operators when acting on the magnetic oscillator functions 

= e~^^^(/)^(x + gA^), where the magnetic length = \J~^- If we neglect the second and 
third terms in the effective Hamiltonian (1), we can describe the low-energy Landau levels 
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as a sequence of levels EEHSO] 



|n"): 




n = 0, 1, 




n>2, 



(12) 



with energies e^c = afiujc\^n{n — 1), where c^c = eS/m, and a distinguishes between the 
conduction (a = 1) and the valence (a = — 1) band. To incorporate the v-^ and the strain- 
induced terms into this picture, we foUow the approach employed before to describe the effect 
of the 73 coupling on the LL structure of graphite ^3Jj. We describe the new eigenstates, 
each of them a linear combination of an infinite number of functions ipn^ with vectors which 
(2n — l)th and (2n)th component create a minimal subspace required to describe the nth 
LL in the absence of v^. In the corresponding infinite Hamiltonian matrix, 
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(13) 



where 



H(n) 



Din) 



'%fvjjc\/n{n + 1) 
7i 




-ii^\/2heBn 



W 




(14) 
(15) 
(16) 



the first term in the low-energy Hamiltonian results in decoupled diagonal 2x2 blocks, the 
second term results in some off-diagonal couplings between functions %l)n and '^nis, while 
the strain-induced third term leads to off-diagonal couplings between functions ipn and il)n±2 
[32] . We then truncate the infinite matrix, restricting the calculation to a given n LLs and 
diagonalise the resulting Hamiltonian numerically. The number of the basis vectors required 
in the calculation in order to describe correctly the low-energy LL structure increases with 
decreasing magnetic field (refiecting growing importance of the '^3 terms at weaker fields). 



11 



0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 



magnetic field [T] magnetic field [T] 



= 








— v=+4 




— v=+4 ^00^ 


— v=+4 >^ 


v=+8 J 






V=+S y/^'V=+12 . 


— v=+i2 






y/^ V=+16 


v=+16 y/^'^^ 


/A 








/ — v=+12 

1 ^ V=+16 ^yf 
1 ✓ 






/a- — v=+12- 

fy^^^ v=+16 
/ . . . . 







^ .... 1^ .... 1 1 . . . . 

0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



magnetic field [T] magnetic field [T] 

FIG. 5: Top row: fan plots of Landau levels for (strained) bilayer graphene, for representative 
points in the ^w) space. Numbers in boxes correspond to the additional degeneracy of a given 
Landau level (valley and spin included) and v on some of the graphs denotes the most stable filling 
factor. Bottom row: activation energies for the quantum Hall effect in strained bilayer graphene 
with various integer filling factors. For a large enough strain, filling factor v — ±4 would be the only 
persistent feature in the low- field quantum Hall effect. The appearance of a local minimum in the 
electron dispersion upon a collision of two Dirac points {w — — 5e*) is manifested by an intermediate 
saturation of the magnitude of the activation gap within the interval 0.2T < B < 0.4T. 

In this work, matrices of the dimension 1200 have been used to produce the spectra, which 
are shown in the top row of Fig. ^foi w = 0^ w = 5e*, w —56* and w — i5e*. 

The additional degeneracy of the 0th LL in bilayer graphene at weak fields is connected 
to the number of Dirac cones in the electronic spectrum. At magnetic fields such that the 
inverse of the magnetic length is smaller than the distance in the reciprocal space between 
a pair of cones, each of them provides additional four- fold degeneracy to the LL (just like 
in monolayer graphene [2l^ l3j). At stronger fields, contributions of separate cones can no 
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longer be resolved, leading to an 8- fold degenerate Oth LL [F^ ^6]. Hence, for the strain w 
within the region shown in yellow in Fig. [3| the zeroth LL becomes 16-fold degenerate at 
very weak but nonzero fields as three LLs, 8-, 4- and 4- fold degenerate, merge together, 
refiecting the existence of four Dirac cones in the spectrum. However, for regions of greater 
strain, shown in Fig |3] with dark and light shading, the zeroth LL is 8-fold degenerate at 
all fields due to only two cones in the electronic spectrum. As the size of the gap between 
two LLs decreases with the LL index, the filling factor u determining the biggest activation 
gap at low fields is closely connected to the additional degeneracy of the zeroth LL. Graphs 
showing the magnitude of activation gaps for filling factors u = ±4, u = ±8, u = ±12 and 
= ±16 as a function of magnetic field for the values of w representing all regions from Fig. 
[3| are shown in the bottom row in Fig. [5| For the region of weak strain, here portrayed 
with w = 0^ the biggest activation gap occurs at = ±8, whereas for the two other regions 
{w = be"" ^ w = —56* and w = — i5e*), filling factor u = ±4 is the most stable. For the region 
of intermediate strain {w = —56*), where a local parabolic minimum exists in the electronic 
dispersion, the activation energy between the Oth and the first LL experiences a very unusual 
intermediate saturation, indicating that one of the LLs gets stuck in that minimum. 

VI. ACKNOWLEDGEMENTS 

This project has been funded by the EPSRC grants EP/G041954 and Science & Innova- 
tion Award EP/G035954. This review is based upon the lecture course "Electronic properties 
of graphene" taught by V. Fal'ko at the 5th Windsor Summer School, 9 — 21/8/2010, in 
Windsor, United Kingdom [33j. 



[1] A.K. Geim, Science 324, 1530 (2009). 

[2] K.S. Novoselov et al, Nature 438, 197 (2005). 

[3] Y. Zhang, Y.-W. Tan, H.L. Stormer, and P. Kim, Nature 438, 201 (2005). 

[4] D.S.L. Abergel, V. Apalkov, J. Berashevich, K. Ziegler, and T. Chakraborty, Adv. Phys. 59, 
261 (2010). 

[5] E. McCann and V.I. Farko, Phys. Rev. Lett. 96, 086805 (2006). 

[6] K.S. Novosleov et al, Nature Phys. 2, 177 (2006). 

13 



[7] I.M. Lifshitz, Sov. Phys. JETP 11, 1130 (1960). 
[8] J.L. Manes, Phys. Rev. B 76, 045430 (2007). 

[9] V.M. Pereira, A.H. Castro Neto, and N.M.R. Peres, Phys. Rev. B 80, 045401 (2009). 
[10] J.D. Bernal, Proc. R. Soc. A 106, 749 (1924). 

[11] X. Du, I. Skachko, A. Barker, and E.Y. Andrei, Nature Nanotech. 3, 491 (2008). 

[12] K.I. Bolotin et al, Sol. St. Commun. 146, 351 (2008). 

[13] B.E. Feldman, J. Martin, and A. Yacoby, Nature Phys. 5, 889 (2009). 

[14] R.T. Weitz, M.T. Allen, B.E. Feldman, and A. Yacoby, Science 330, 812 (2010). 

[15] J. Martin, B.E Feldman, R.T. Weitz, M.T. Allen, and A. Yacoby, Phys. Rev. Lett. , (2010). 

[16] E.V. Castro, N.M.R. Peres, T. Stauber, N.P. Silva, Phys. Rev. Lett. 100, 186803 (2008). 

[17] F. Zhang, H. Min, M. Polini, and A.H. MacDonald, Phys. Rev. B 81, 041402 (2010). 

[18] O. Vafek and K. Yang, Phys. Rev. B 81, 041401 (2010). 

[19] Y. Lemonik, I.L. Aleiner, C. Toke, and V.I. Fal'ko Phys. Rev. B 82, 201408 (2010). 

[20] R. Nandkishore and L. Levitov, Phys. Rev. Lett. 104, 156803 (2010). 

[21] V.N. Kotov, B. Uchoa, V.M. Pereira, A.H. Castro Neto, and F. Guinea, arXiv.:10123484. 

[22] Y.-W. Son, S.-M. Choi, Y. P. Hong, S. Woo, S.-H. Jhi, |arXiv:1012.0643 (2010). 

[23] J.C. Slonczewski and P.R. Weiss, Phys. Rev. 109, 272 (1958). 

[24] M.S. Dresselhaus and G. Dresselhaus, Adv. Phys. 30, 139 (1981). 

[25] J.R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966). 

[26] D.S.L. Abergel and V.I. Fal'ko, Phys. Rev. B 75, 155430 (2007). 

[27] T. Ando, J. Phys. Soc. Jpn. 76, 104711 (2007). 

[28] M. Koshino and T. Ando, Phys. Rev. B 77, 115313 (2008). 

[29] M. Nakamura, L. Hirasawa, and K.-I. Imura, Phys. Rev. B 78, 033403 (2008). 

[30] M. Mucha-Kruczynski, D.S.L. Abergel, E. McCann, and V.I. Fal'ko, J. Phys. Condens. Mat. 

21, 344206 (2009). 
[31] K. Nakao, J. Phys. Soc. Jap. 40, 761 (1976). 

[32] Note, that two rows and columns in the above matrix, identically equal to zero, give rise to 
two solutions at zero energy, which correspond to unphysical eigenstates (zero vectors) and 
should not be confused with the zero energy Landau levels. 

[33] V.I. Fal'ko, Lecture notes on "Electronic properties of graphene" , 5th Windsor Summer School, 
9 - 21/8/2010, Windsor, United Kingdom. 

14 



